Astron. Nachr. / AN 999, No. 88, 789-|792l(2010) / DOI please setDOI! 



Asteroseismic study of 7 Doradus members of the cluster NGC 2506 

A. Grigahcene^ *, A. Moya^, J.-C. Suarez'^ and S. Martin-Ruiz^ 

^ Centro de Astroffsica, Faculdade de Ciencias, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

^ Laboratorio de Astroffsica Estelar y Exoplanetas, LAEX-CAB (INTA-CSIC), PO BOX 78, 28691 Villanueva de la 

Canada, Madrid, Spain 
^ Instituto de Astroffsica de Andalucfa (CSIC), CP3004, Granada, Spain 



O 

(N 

Ph. 
< 



6 



> 

O 

o 



Received 2010, accepted 2010 
Published online later 

Key words Editorial notes - instruction for authors 

The present work performs an asteroseismic study of a 7 Doradus star member of the open cluster NGC 2506. This works 
conjugates the observational information of the star and the strong constraints imposed by the cluster memebership, with 
the latest developed theoretical asteroseismic tools for the modelling of this type of objects, in particular the Frequency 
Ratio Method (FRM) and Time-Dependent Convection (TDC), to study the 7 Dor stars found in the NGC 2506 stellar 
cluster. The use of both techniques gives us the opportunity of constructing a self-consistent procedure, allowing the mode 
identification and improving the modelling of the 7 Dor members of the cluster. In this work we present the result of the 
analysis of the first target of our project, showing the advantadge of modelling cluster member's 7 Dor stars. This is the 
first step in a more ambitious project. In the near future, we plan to expand the study to other variable members of the 
NGC 2506 cluster. This will improve our knowledge of 7 Dor pulsatorsand probe the self-consistency of our procedures. 

© 2010 WILEY- VCH VcrlagGmbH&Co.KGaA, Wcinhcim 



1 Introduction 

The 7 Doradus (Dor) stars are suitable objects for astroseis- 
mic study. They are high-order gravity (g) mode pulsators, 
which permits the probing of their stellar coreand the use 
of the asymptotic approximation, for w hich analytical so 



lution s of the oscillation equations exist (ISmevers & Mova 



2005 



2007 ) The Frequency R atio Method (FRM; 



Mova et al 



Suarez et al.ll2005l) . based on this property, is partic- 



ularly adapted for obtaining asteroseismic information on 
7 Dor pulsating stars that show, at least, three pulsation fre- 
quencies. The method provides an estimate of the identifi- 
cation of the radial order n and the spherical degree £ of 
observed frequenciesand of the integral of the buoyancy 
frequency (Brunt-Vaisala) weighted over the stellar radius 
along the radiative zone (/obs)- 

The excitation mechanism of 7 D or stars was a n out- 



stan ding problem until the works of iGuzik et al. I (l2000h 



and iDupret et al.l (120051) who used respectively, frozen and 
Time-Dependent Convection (TDC) theories showed that 
the position of the base of the convective envelope is the 
key for driving 7 Dor g modes. Moreover, TDC is formu- 
lated in the framework of mixing length theory and there- 
fore the 7 Dor instability strip is predicted to be sensitive to 
the mixing-length parameter auur- 

The asteroseismic modelling of 7 Dor, even armed with 
these techniques, still needs the global parameters of the 
observed stars. One of the ways to obtain a significant en- 
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hancement is to use the strong constraints imposed by the 
cluster membership for examp le, the ag e, metallicity and 
distance dCreevev et al. 2009; Suarez et al. 2002, 2007). 

At the beginning of the recognition of the 7 Dor stars as 
a new class of variable stars, a great effort has been devoted 
to their detection in clusters. At that time the main aim of the 
community was to check whether this type of pulsations was 
linked or not to age and obviously, ste llar clusters were the 
most suitable targets for this purpose (IMartm & Rodriguez 
2002h . 

One of the richest cluster in 7 Dor members is the 
old open cluster NGC 2506, for which photometric ob- 
servations provide a considerable number of variable stars 
dArentoft et al.ll2007h . Among them, 15 have been identi- 
fied as 7 Dor stars candidates (labeled as Vn - V25 in that 
paper). 

This work attempts to investigate the impact of using 
a procedure which includes FRM and TDC on the mod- 
elling of 7 Dor, stars taking advantage of their membership 
in clusters. In particular, we here provide preliminary results 
for the 7 Dor star Vn, belonging to the open cluster NGC 
2506. 



2 Observational Data 

The old open cluster NGC 2506 (a, d)2ooo=(08h OO'" OP, - 
10° 46' 12") has been reported to hold around 20 oscillating 
stars of different types, 15 of them being 7 Dor candidates. 
The physical properties of the cluster are: distance = 3460 
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Fig. 1 HR Diagram position of 7 D or variables. Crosse s 
represent the stars parameters giv en by Henry et^(2007 ). 
The NGC 2506 7 Dor members dArentoft et alj|2007l) are 
displayed by filled circles, among which Vn (labelled with 
'11'). Dashed lines show the the oretical 7 Dor insta bility 
strip for q;mlt=2 as calculated by iDupret et al.l (l2005b . 



Obviously, stellar models could differ considering the spe- 
cific values chosen for these global parameters (Toff, Metal- 
licity. Log g, etc.) or for the values of the stellar modelling 
such as the Mixing-Length parameter aMLT or overshoot- 
ing Qfov, etc. 

We have constructe d a grid of models, using the CE- 



SAM equilibrium code (iMorel & Lebreto 



20081), varying 
steps of 
0.12, 



the mass (in the range [1.25, 2.10] M© with 
O.OIMq), the metallicity (with values [M/H]=0.08, 
— 0.32and —0.52), the mixing-Length parameter (values 
«MLT= 0.5, land 1.5) and overshooting (values aov= 0.1, 
0.2and 0.3). Standar d physical inputs for 7 Do r and S Scuti 



stars are used (see ICasas et al.l 120061 120091 for details) 



The pulsational observabl es have been calculated using 
the GraCo pulsational c ode ( Mova. Garrido & Dupretl2004 ; 
Mova & GairidolbOOSb . Both codes (CESAM and GraCo) 
have been part of the ESTA group for c omparison of codes 



have been part or the fca lA group tor c ( 
(iLebreton et al.ll2008t iMova et al.KOOSh 



The number of models lying in the error box and ful- 
filling the cluster membership (sub-solar metallicity and 
age > 1800 Myr), is huge (1210 in our grid). The mini- 
mum and maximum values obtained are displayed in Ta- 
ble [T] The sub-solar metallicity constraint is automatically 
fulfilled when the age restriction is imposed, that is, there 
are no models of solar metallicity with an age larger that 
1800 Myr and the global parameters of NGC 2506 Vn. This 
agrees with the global determination for the cluster 



Table 1 Minimum and maximum values of the physical 
parameters of NGC 2506 Vn matching its observational pa- 
rameters, with ages larger than 1800 Myr 



M 


Tcff 


[Fe/H] 


L 


Age 


Qmlt 




(Mq) 


(K) 


(dex) 


(L0) 


(Gyr) 






1.25 


6930 


-0.52 


6.03 


1.80 


0.5 


0.1 


L45 


7229 


-0.12 


6.92 


3.12 


1.5 


0.3 



PC age > 2 Gyrand a metall icity of [Fe/H] = -0.20 ± 0.1 
dex (see lArentoft et al. (2007) and references therin). 

In Fig. [T]we show these 7 Dor candidates of NGC 2506 
(filled circles) together with known bonafide (crosses) 7 Dor 
stars in the HR diagra m. We also show the theoretical insta- 
bility strip obtained by lDupret et al.l (120051) . 

In this study we concentrate on NGC 2506 Vn (V = 
15.454) which is labeled '11' in Fig.[T] Its global parameters 
are given by Tcff = 7080 ± 150 K and Log (L/Lq) = 0.81 
± 0.03. It has been chosen as a first target because it shows 
the highest frequencies: 1.165, 1.270 and 1.400 (d"^), re- 
spectively. 



3 Modelling of NGC 2506 Vn using global 
parameters 

The first step in the modelling of Vn was the selection of 
those models matching the global parameters of the star. 



4 Theoretical Results - FRIVI 

The next step in our study was the inclusion of the astero- 
seismic constraints. To do so, we followed a self-consistent 
procedure presented by Moya et al. (2008) used in different 



asteros ei smic studies of 7Doradu s stars (iRodrfguez et al 
2006aB lUytterhoeven et al.ll2008h and recently applied to 
the determination of the physical characteristics of the 
first planetary syste m observed by direct imaging, HR 8799 
dMova et al.ll2010allbl) . 

We first applied the FRM. This method provides sets 
of possible mode identifications (n, £) for the three ob- 
served frequencies and the corresponding value for the 
Brunt-Vaisala integral (lobs). Adding this new constraint, 
the number of accepted models is reduced by 70 %. Then, 
we searched for those models matching the observed fre- 
quencies with the predicted mode identification and Jobs- 
The model selection was done using statistics. The total 
number of modes selected is 50, for which the standard de- 
viation of the frequency fit is lower than 0.5 yuHz. Table |2] 
shows the models with the best ^hd the mode identifica- 
tion predicted. 

5 Theoretical Results - TDC 

The second step of our procedure is the study of the en- 
ergy balance of the theoretical modes using TDC. We have 
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Table 2 Best models fitting the observed frequencies, with the mode identification and Jobs predicted by the FRM. 



M/Mq 


Model 
[Fe/H] 






ni 


^1 


n2 


^2 


na 






1.38 


-0.32 


1.5 


0.1 


-50 


2 


-46 


2 


-24 


1 


1.598E-06 


1.36 


-0.32 


1.5 


0.1 


-49 


2 


-45 


2 


-41 


2 


2.563E-06 


1.27 


-0.52 


1.5 


0.2 


-30 


1 


-48 


2 


-25 


1 


6.073E-06 


1.36 


-0.32 


1.5 


0.2 


-48 


2 


-44 


2 


-23 


1 


6.861E-06 


1.28 


-0.52 


1.5 


0.1 


-56 


2 


-30 


1 


-47 


2 


8.374E-06 


1.33 


-0.32 


1.5 


0.2 


-48 


2 


-44 


2 


-40 


2 


9.794E-06 



[Fe/H]=-0.52 t 
[Fe/H]=-0.32 □ 



1.38 

1.36 

1.34 
a 1.32 
S 1.3 

1.28 

1.26 

1.24 

1800 2000 2200 2400 2600 2800 3000 3200 
Age in Myr 

Fig. 2 Age - Mass diagram of the models fulfilling all the 
observational constraints. The metallicity of each model is 
also shown. 

accepted those models selected in the previous step that pre- 
dict over-stable observed frequencies. The models fulfilling 
all the asteroseismic constraints are drastically reduced to 
those displayed in Fig. |2] In this figure we show the Age - 
Mass values of these models together with their metallicity. 
The minimum and maximum values obtained are displayed 
in Table [3] 

A first consequence of this study is that models with 
an age around 1800 - 2000 Myr have an internal metallic- 
ity similar to that observed for the cluster ([Fe/H]=-0.32), 
masses in the range [1.33,1.38] Mq and radius around 1.77 
Rq. Lower metallicity models can also fulfill all the obser- 
vational constraints but with ages larger than expected for 
the stars ofNGC 2506. 

6 Conclusions 

With this prospective study we showed that the application 
of FRM and TDC to the cluster member 7 Dor candidate 
NGC 2506 Vii enhances significantly our determination of 
the mass, radius and age. Specifically, its metallicity is con- 
firmed to be sub-solar. The energy transport efficiency of the 
outer convective zone is high (omlt = 1-5). These results 



must be confirmed with the extension of this study to other 
7 Dor stars of this cluster. 

This is a first step of a more ambitious project. In the 
near future, we plan to expand the study to other variable 
members of the NGC 2506 cluster. A better knowledge of 
variable membership jOttosenl l2008b and their global pa- 
rameters will improve our knowledge on 7 Dor pulsators 
and probe the self-consistency of our procedures. 
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Table 3 Minimum and maximum values of the physical parameters of NGC 2506 Vn fulfilling all the constraints 



Mass 




Teff(K) 


[Fe/H] 


L/L0 


R/R0 


Age (Myr) 


omlt 




1.25- 


1.29 


6997 


-0.52 


6.07 


1.6 


1832 


1.5 


0.1 


1.33- 


1.38 


7223 


-0.32 


6.83 


1.774 


2585 


1.5 


0.3 
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